Single Transverse-Spin Asymmetry in Very Forward and Very Backward Neutral 
Particle Production for Polarized Proton Collisions at = 200 GeV 
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In the 2001-2002 running period of the Relativistic Heavy Ion Collider (RHIC), transversely 
polarized protons were accelerated to 100 GeV for the first time, with collisions at y/s = 200 GeV. 
We present results from this run for single transverse spin asymmetries for inclusive production of 
neutral pions, photons and neutrons of the energy region 20 - 100 GeV for forward and backward 
production for angles between 0.3 mrad and 2.2 mrad with respect to the polarized proton direction. 
An asymmetry of ^jv = (-0.090 ± 0.006 ± 0.009) x (l.Ol^I) was observed for forward neutron 
production, where the errors are statistical and systematic, and the scale error is from the beam 
polarization uncertainty. The forward photon and 7r°, and backward neutron, photon, and ir° 
asymmetries were consistent with zero. The large neutron asymmetry indicates a strong interference 
between a spin-flip amplitude, such as one pion exchange which dominates lower energy neutron 
production, and remaining spin non-flip amplitudes such as Reggeon exchange. 

PACS numbers: Valid PACS appear here 



Polarized proton beams were accelerated to 100 GeV 
for the first time in the 2001-2002 running period of the 
Relativistic Heavy Ion Collider (RHIC). We present re- 
sults from an experiment at the 12 o'clock intersection 
point (IP12), for the single-spin asymmetry An for colli- 
sions of transversely polarized protons with unpolarizcd 
protons, producing very forward (and very backward) 
neutral particles. The goal of the measurement was to 
search for a production process with a significant asym- 
metry that could be used as a local polarimeter at the 
experiments, to monitor the polarization direction in col- 
lision at RHIC. The very forward region was chosen due 
to accessibility with the existing RHIC detectors. The 
experiment was designed to search for photon and ir° 
asymmetries, based on the observation of a large pion 
asymmetry for yfs = 19.4 GeV.[l[ However, a large neu- 
tron asymmetry was observed, which is reported here. 

The measured cross sections for inclusive production of 
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neutrons at large x f for unpolarized proton-proton colli- 
sions are consistent with one pion exchange (OPE) model 
predictions from = 7 to 64 GeV.@,Hll,lS] The model 
has been used to describe exclusive diffractive produc- 
tion of neutrons [f| , production near the exclusive limit 
at large xf, and inclusive neutron and proton production 
away from this limit (see 0). More recently, with inter- 
est in HERA e-p collider studies of proton fragmentation, 
OPE is used to describe proton and photon induced pro- 
duction of neutrons which may probe the pion structure 
function at small x.{5, 7} For polarized phenomena, Soffer 
and Tornqvist Q applied the OPE model to describe the 
polarization of inclusive A production in p-p collisions. Q 
A single transverse spin asymmetry, the analyzing power 
An, must arise from an interference of spin flip and spin 
non-flip amplitudes. The OPE amplitude is fully spin 
flip; a spin non-flip amplitude contribution can be de- 
scribed in Regge theory as due to Reggeon and Pomeron 
exchange. [l(| To the extent that the cross section is de- 
scribed well by spin-flip OPE, An is a sensitive test of 
this picture, and is sensitive to a small admixture of a 
non-spin flip contribution at the amplitude level. 

The analyzing power is the azimuthal modulation of 
the cross section, or left-right asymmetry of a vertically 
polarized beam, of produced particles relative to the spin 
direction of a transversely polarized beam on an unpo- 
larized target. The measurement consists of observing 
the azimuthal dependence of the production of identified 
particles, normalized by an independent measurement of 
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the beam polarization. 

The RHIC polarized proton beams were vertically po- 
larized, and bunched, with a bunch length of ~3 ns, with 
the bunches separated by 212 ns. The 55 bunches in one 
ring (the Blue ring) had alternating spin directions in suc- 
cessive bunches, and the other ring (the Yellow ring) had 
alternating spin directions for successive pairs of bunches. 

RHIC timing signals identified the bunches and their 
polarization sign. This provided collisions between pro- 
tons with all combinations of vertical spin orientation, at 
the same time. Asymmetries could be measured for ei- 
ther beam polarization, corresponding to very forward or 
very backward production, by averaging over the polar- 
ization directions of the opposing beam. After averaging 
over the spin directions for one beam, and over many fills 
of RHIC, the residual beam polarization was less than 1% 
of the RHIC beam polarization. A typical fill of RHIC 
for this first run had 3 x 10 10 protons in each bunch. 
With (3* — 3 m at the IP 12 collision location, the typical 
luminosity was L = 1 x 10 30 cm _2 s _1 and the average 
beam divergence was 0.1 mrad rms. 

The polarization of each beam was measured by scat- 
tering from ultra-thin carbon targets that were inserted 
into the beams for the measurements, at a location 
where the beams are separated. 11] Measurements took 
less than one minute, giving a polarization sensitivity of 
about 3%, and were typically taken every two hours of a 
six hour store. The measured beam polarizations for the 
2001-2002 run were P Blue = (0.11±0.002±0.02) x (1.0± 
0.32), and Pyellow = (0-16 ± 0.002 ± 0.02) x (1.0 ±0.32). 
The first uncertainty is statistical, the second the sys- 
tematic on the asymmetry measurement, and the scale 
error is from the calibration of the polarimeter analyzing 
power. [H [3 

The experimental apparatus at IP 12 consisted of a pair 
of scintillator hodoscopes placed forward and backward 
of the collision point that identified proton-proton col- 
lisions, and detectors placed in the very forward direc- 
tions, centered on a production angle of 0° downstream 
of the RHIC DX dipole magnets, which removed pro- 
duced charged particles (Fig. [T|). The collision-trigger 
hodoscopes (beam-beam counters) were located ±1.85 m 
from the IP center. The hodoscopes were formed with 
four sets of rectangular scintillators, with full azimuthal 
coverage, and pseudorapidity acceptance | ^7 1 = 2.2 to 3.9 
in the vertical and horizontal directions from the beams. 
The time resolution of the hodoscopes gave a vertex res- 
olution of 23 cm. The vertex distribution had a 54 cm 
rms, and non-collision events within the collision trigger 
were estimated to be 3%. 

The east detector consisted of a scintillation counter 
charged particle veto (dimensions 10.5 cm x 25 cm x 
0.6 cm, CV in Fig. []}, a preshower electromagnetic 
calorimeter, and a 60-element array of PbWC>4 crystal 
electromagnetic calorimeter (EMCal). These were fol- 
lowed by scintillation counters on either side of an ab- 
sorber to observe shower leakage. The preshower and 
EMCal arrays were built from PbW04 crystals 2.0 x 
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FIG. 1: A plan view of the experiment at IP12, not to scale. 
Shown are the principal components of the experiment: the 
beam-beam collision trigger hodoscopes, the electromagnetic 
calorimeter (east) facing the Blue beam, and the hadronic 
calorimeter (west) facing the Yellow beam. 



2. Ox 20.0 cm 3 , coupled to 3/4-inch photomultiplier tubes 
(HPK R4125) via silicon cookies through filters of 1/10 
light reduction. The preshower consisted of five crystals 
forming a horizontal hodoscope, 10 cm (horizontal) x 20 
cm (vertical) x 2 cm (deep). The EMCal was a 5 (hor- 
izontal) x 12 (vertical) crystal array, 20 cm (22 radia- 
tion lengths) deep. The preshower and EMCal were cali- 
brated in the electron beam of the Final Focus Test Beam 
Facility at SLAC. The energy resolution was AE/E ~ 
10%/a/-E(GcV) and the position resolution was 0.1 cm, 
independent of the hit position. The energy resolution 
was confirmed by the width of the observed 7r° peak dur- 
ing the experiment. For hadronic showers, the position 
resolution was estimated to be 0.5 cm using the simu- 
lation program Geant3[14], which had been tuned to 
match the EMCal response to the electrons. However, 
the EMCal was 1 interaction length with an estimated de- 
posited energy fraction for neutrons of Ed ep ■ E n ~ 1 : 3 
with a > 100% energy resolution. The EMCal was fol- 
lowed by two scintillators that covered the area of the 
EMCal (dimensions 10 cm x 24 cm x 0.6 cm), with a 
2.8 cm thick iron block in between, Nl and N2 in Fig. Q] 
Data were taken with the primary trigger requirement 
of a collision, with > 1 minimum ionizing particles in each 
beam-beam hodoscope, within a 5 ns collision window, 
and with a minimum energy deposited in the EMCal, 
>5 GeV. Analog to digital converters recorded the ener- 
gies for each calorimeter block, time to digital converters 
recorded the summed signals from the beam-beam coun- 
ters, and each counter was latched. The polarization di- 
rection information (up or down) was provided by the 
RHIC accelerator, and the recorded data were sorted by 
the polarization directions of the Blue and Yellow beams. 
The trigger rate was 100-300 Hz, and the average live 
time of the experiment was 75%. 70 million events were 
collected. 

Events were reconstructed in the EMCal by combining 
clusters of hits to provide a total observed energy, and 
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with positions of the showers obtained from the energy- 
weighted average of the positions of the hit blocks. 

Photon samples and neutron samples were identi- 
fied using the observed preshower and postshower en- 
ergy. The photon samples were identified as events 
with no charge veto hit, preshower energy higher than 
15 MeV, and no postshower activity ("7-ID"). The 
neutron-sample identification required no charge veto hit, 
no preshower energy, and postshower activity ("n-ID"). 
Both are triggered by the collision requirement and en- 
ergy in the EMCal, >5 GeV. Fig.[5]shows the distribution 
of energy deposited in the EMCal for the photon-sample 
(a) and neutron-sample (b) selection requirements. Note 
that for Fig. [5] (b) , the spectrum shape is predominantly 
the result of the poor energy resolution of the EMCal for 
neutrons. 

The fraction of photons and neutrons in each sample 
was obtained by comparing the fractions of events in the 
samples with and without preshower activity. These frac- 
tions can be calculated based on the number of real pho- 
tons and neutrons in each sample, using the preshower 
conversion probabilities for photons and neutrons. For 
the fraction of neutron conversions in the preshower that 
have post shower activity, to the fraction that convert 
in the EMCal with observed energy in the post-shower, 
we have assigned the relative probability of 0.5 ± 0.5. 
The fraction of neutrons with energy >20 GeV in the 
7-ID sample is 0.18 ± 0.07 (the purity of photons is 
P 7 = 0.82 ± 0.07); the fraction of photons in the n- 
ID sample is 0.006 ± 0.005 (the purity of neutrons is 
P n = 0.994 ± 0.005). 

Fig. [3] shows the invariant mass for events with two 
clusters of energy in the EMCal, for no additional require- 
ment for photon identification, and for the 7-ID sample. 
The 7T° peak mass shifts lower by 5 MeV when no pho- 
ton identification is required due to energy lost in photon 
conversion upstream; both the mass shift and reduction 
in number of events are understood from simulation us- 
ing the upstream geometry, in particular conversion in 
a 2.1 cm thick stainless steel plate just downstream of 
the DX magnet shown in Fig. [TJ The it peak center 
for events with no photon identification was used for an 
absolute energy calibration and tower-by-tower final cal- 
ibration. To measure the analyzing power for ir° events, 
we use the larger data sample without photon identifica- 
tion, selecting events within ±20 MeV of the peak center, 
with 21% ± 3.8% background, estimating the background 
using a Gaussian + polynomial fit. 

To measure a left-right asymmetry, it is necessary 
to eliminate 0° production. We required the observed 
shower center to have radii 5-40 mm from the center of 
the detector, or production angles 0.3 - 2.2 mrad. The 
EMCal hits were assigned coordinates (r,<^>) about the 0° 
reference axis, and binned in (f> counterclockwise, with <p 
= along the vertical axis. A vertical beam polarization 
would give a left-right asymmetry for <f> — 7r/2, and a 
radial beam polarization would give an asymmetry for <j> 
= and <fi = it. 
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FIG. 2: Panel (a) shows the deposited energy distribution 
in the EMCal for events used to select the photon sample, 
for events with no charged veto, and no postshower energy. 
The solid curve displays all events, the dashed curve dis- 
plays events with preshower energy (this is the 7-ID sample), 
and the dotted curve shows events with no preshower energy. 
Panel (b) shows the energy distribution for the events used 
to select neutrons, for events with no charged veto, and with 
postshower activity. The solid curve is for all events, the 
dashed curve is for events with no preshower energy (n-ID 
sample), and the dotted curve is for events with preshower 
energy. 




FIG. 3: Invariant mass of pairs of energy clusters in the 
EMCal, for 444K events with no additional selection require- 
ments (solid) and for 35K events with photon identification 
(dashed). 
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were calculated by: 



FIG. 4: Azimuthal dependence of asymmetry for the n-ID 
sample produced forward with respect to the polarized pro- 
ton direction, based on the east detector. The error bars are 
statistical. 



To calculate the asymmetry, we used a method of 
analysis that reduces systematic errors from detector 
and luminosity asymmetries. We defined the geomet- 
ric mean for the number of observed events with beam 
proton spin up, producing particles to azimuthal an- 
gle 0, -/Vf^, and the number of events with beam spin 
down, producing particles to azimuth (f) + ir, Ni^+k, 
U Ni^Nl^+k. A corresponding geometric mean was de- 
fined for events with spin down, to angle 0, and events 
with spin up, to angle + tt. The asymmetry is the dif- 
ference of the geometric means, divided by the sum, and 
normalized by the measured beam polarization Pb'. 



A((j)) = 1 y/N^NL^-y/Nij+yNM 

We define An as the amplitude of the asymmetry mod- 
ulation: 



A(4>) = A N sin( 



>), 



(2) 



where 0o allows a deviation of the polarization direction 
from = 0. Fig. [4] shows the azimuthal dependence of 
asymmetry for the n-ID sample for forward production 
with the Blue beam polarized and summing over the spin 
states of the Yellow beam. A significant non-zero asym- 
metry A^f was observed with 0o = —0.15 ±0.07 and a 

reduced \ 2 = 1-5- Similarly the asymmetries A N ~^ for 
the n-ID sample for backward production (Yellow beam 
polarized, averaging over spins of the Blue beam) and 

A N for the 7-ID sample for forward and backward 
production were obtained. Polarized Blue beam results 
used 0o = —0.15, and polarized Yellow beam results used 
00 =0. 

The asymmetries of neutrons (A N ) and photons (A N ) 
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where P n and P 7 are the purity of neutrons in the n-ID 
sample and that of photons in the 7-ID sample respec- 
tively. 

To confirm this neutron asymmetry signal during the 
run, we added a partial hadronic calorimeter (HCal) fac- 
ing the Yellow beam, the west detector in Fig. [TJ The 
setup included a photon veto (7 in Fig. [I}, consisting 
of a 5-cm thick lead block followed by a scintillator, 
the hadron calorimeter, and five PbWC>4 crystals that 
were set up as a horizontal hodoscope to identify left 
vs. right production. The HCal had a transverse dimen- 
sion of 10 cm x 10 cm, and a length of 2 interaction 
lengths. The PbWC>4 crystals were the same dimension 
as used for the east detector. The position resolution of 
the shower center was estimated to be 3 to 4 cm, from 
a Geant3 simulation. [l4T| The asymmetry was obtained 
by comparing the left and right scattering only (the az- 
imuthal dependence was not available with this setup), 
and was corrected by the factor 0.69 for the integration 
over the azimuth, due to the azimuthal dependence of 
the analyzing power for vertical beam polarization. For 
this measurement, neutrons with a shower maximum in 
the central crystal were not included, and a HCal en- 
ergy of 28 GeV or higher in the trigger and 30 GeV 
in the offline analysis was required. The energy reso- 
lution of the HCal was 40-50% for energy greater than 
20 GeV. Thus, the neutron event selection was matched 
as far as possible for the HCal and EMCal setups. The 
measured forward asymmetry with the Yellow beam po- 
larized, and summing over the spin states of the Blue 
beam, was An — —0.135 ± 0.018, a significant and large 
asymmetry, and consistent with the result from the east 
detector at the 2-tr level after including systematic un- 
certainties. 

The HCal was a prototype of the Zero Degree 
Calorimeters (ZDC) used at RHIC, one third of the 
length of the RHIC ZDCs.0 With only two interac- 
tion lengths, the measured energy can range from a small 
fraction of the actual neutron energy for conversions near 
the downstream end of the HCal, to the full neutron en- 
ergy. The neutron asymmetry reported here has been 
used in subsequent RHIC runs by the PHENIX exper- 
iment to set the proton spin direction at collision. A 
neutron energy spectrum has been reported [T3|, which 
is consistent with the ISR results Q , in particular show- 
ing a peak in the energy spectrum at about Xf = 0.8. 
For the experiment reported here, the range in px for 
the neutrons can only be estimated from the limits of 
the acceptance, and the minimum energy requirement, 
which give 0.01 GeV/c < p T < 0.22 GeV/c, for the lim- 
its 30 GeV < E n < 100 GeV and 0.3 mrad < 6 n < 2.2 
mrad. Greatly improved information on the kinematical 
variables will be forthcoming from the PHENIX experi- 
ment. 
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For both calorimeters, the identified neutron sample 
may include K\. At ISR energies, the K\ fraction to 
neutrons for a similar kinematical region was estimated 
to be 3-4%, from observed charged kaon samples. 3 We 
have included no correction for K\ background in the 
neutron results. 

Table Q] gives the corrected asymmetries from the east 
detector for neutron, photon, and 7r° forward produc- 
tion and for backward production. The systematic uncer- 
tainties include contributions from scattered beam back- 
ground, particle misidentification, and from smearing of 
the hit position in the EMCal. Scattered beam back- 
ground refers to an observed higher rate in the EM- 
Cal towers near the beam pipe of the outgoing RHIC 
beam. Whether or not the background has a spin de- 
pendence, it affects the azimuthal dependence of Eq.([T]). 
Its contribution to the systematic uncertainty was esti- 
mated from the reduced x 2 = 1.5 of the azimuthal fit to 
•A n _Try (((>) , where we have assigned a systematic uncer- 
tainty of %/0.5 x 5 An (statistical) . The scattered beam 
background uncertainties for the photon and ir° are neg- 
ligible. The particle identification uncertainties were pre- 
sented earlier. A hit smearing correction of 0.92 ± 0.08 
was evaluated using a Geant3 simulation of the EM- 
Cal which described the test beam result for position 
resolution. Finally, the beam polarization uncertainty, 
presented earlier, scales both the asymmetry and uncer- 
tainties presented in Tabled! preserving the measurement 
significance. The scale factor is (l.O+o;^). 

All asymmetries use events with energies 20 - 100 GeV 
in the EMCal and production angles 0.3 - 2.2 mrad. We 
observe asymmetries consistent with zero for photon and 
7T°, for both forward and backward production, and for 
neutron backward production, and a significant non-zero 
asymmetry for forward neutron production. We have also 
measured the asymmetry for the photon sample obtained 
from conversions near the DX magnet, where the neutron 
contamination is estimated to be less than 12%. This 
asymmetry is 0.016 ± 0.020, consistent with zero. 

In Fig. [5] we present the uncorrected multiplicity ob- 
served by the beam-beam counters for the n-ID sam- 
ple of the EMCal data with deposited energy >20 GeV, 
where the trigger required > 1 hit both forward and 
backward. The vertical and horizontal slats for each ho- 
doscope overlap over 67% of the area, so that a correct 
estimate of the multiplicity is ~ (uncorrected multiplic- 
ity) / 1.67. As seen in the figure, the multiplicity for 
the forward beam counter, in the direction of the neu- 
tron, is low, (multiplicity (forward)) ~ 2. The multiplic- 
ity in the backward beam counter is large, (multiplicity 
(backward)) ~ 7. Therefore, we observe a clear separa- 
tion of beam and target fragmentation multiplicity, with 
a large asymmetry for neutrons produced forward from 
the polarized beam, in the direction of low multiplicity. 
The backward multiplicity is large. This pattern was 
confirmed with the HCal data, where the forward and 
backward beam counters are reversed. The low forward 
multiplicity may be consistent with a simple exchange 
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FIG. 5: Uncorrected multiplicity observed by the beam-beam 
counters for the n-ID sample of the EMCal data with de- 
posited energy >20 GeV. The solid line shows that for the 
forward counter (east) and the dotted line shows that for the 
backward counter (west). 





forward 


backward 


neutron 
photon 

^° 


-0.090 ± 0.006 ± 0.009 
-0.009 ±0.015 ±0.007 
-0.022 ± 0.030 ± 0.002 


0.003 ± 0.004 ± 0.003 
-0.019 ±0.010 ±0.003 
0.007 ±0.021 ±0.001 



TABLE I: Asymmetries measured by the EMCal. The er- 
rors are statistical and systematic, respectively. There is an 
additional scale uncertainty, due to the beam polarization un- 
certainty, of (l.Ot^'al). 



process producing the forward neutron, p-tn, although 
the acceptance of the forward beam counter, with 43 
mrad < 8 < 297 mrad, covers much larger angles than 
the opening angle of, for example, A + — > n ± n + , where 
$open < 3 mrad for an 80 GeV neutron. OPE model de- 
scriptions of neutron production do not in general predict 
the backward multiplicity. 

In summary, we have measured the single transverse- 
spin asymmetry, An, for photon, 7r°, and neutron pro- 
duction at very forward and very backward angles in 
pt+P collisions at y/s — 200 GeV. The asymmetries for 
photon and ir° were consistent with zero, whereas a sig- 
nificant asymmetry was observed for forward neutrons. 
The large neutron asymmetry is accompanied by a small 
forward multiplicity and a large backward multiplicity. 
The asymmetry for neutrons produced backward from 
the polarized beam was consistent with zero. This newly 
discovered asymmetry can be used for a non-destructive 
polarimeter to monitor the beam polarization in colli- 
sion and the physics can be explored with a dedicated 
hadronic calorimeter. 
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